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Abstract 



If} [ The problem of primordial nucleosynthesis of helium in brane cosmology is stud- 

ied using a semi-analytical method, where the Fermi-Dirac statistic is ignored. This 
semi-analytical method agrees with a more complete numerical calculation with a 
C<*) , precision of order of 10% or better. The quadratic term for the matter density 

is the only source considered in the modified Einstein equations predicted by the 
brane configuration. This hypothesis is justified a posteriori. An agreement be- 
tween theoretical and observational values for the helium abundance is obtained if 
the fundamental mass scale in five dimensions is of the order of M ~ 5 x 10 3 GeV . 



PACS numbers: 98.80.Ft, 26.35. +c 



1 Introduction 



The brane model is one of the most interesting cosmological scenario proposed in recent 
times 00 El- For a review of brane cosmology, see references jHEj. In its most simple 
formulation, a five dimensional space-time is considered where a three-brane configuration 
is settled out. The ordinary matter and fields are restricted to "live" on the brane. 
Only gravity and a cosmological constant term are allowed to "live" in the entire four 
dimensional space, called the bulk. One fundamental aspect of this construction is the 
fact that the fifth dimension may be compact but not small: This allows to solve the 
hierarchy problem, explaining why the Planck mass in our usual four dimensional space- 
time takes the huge value Mp\ ~ 10 19 (jeV. In fact, the Planck mass Mpi is related with 
the five dimensional fundamental mass scale M through the relation 

M 2 Pl = M 3 R , (1) 
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where R is the characteristic dimension of the fifth dimension. Hence, the bulk gravita- 
tional coupling can have values compatible with other physical coupling constants pro- 
vided the scale R is very large, inducing at same time a huge value for the Planck mass 
on the brane. In a very fashion variation of this model, Randall and Sundrum [HI Hj had 
considered a five dimensional anti-deSitter model (a negative cosmological constant living 
in the bulk), where gravity has the usual features on the brane even if the fifth dimension 
is infinite. 

One of the distinguishing features of brane cosmology is the presence of a quadratic 
matter term as source of the Einstein equations which describes the gravitational dynam- 
ics on the brane. This quadratic term comes from the matching conditions of the five 
dimensional Einstein equations in presence of a brane, together with the assumption that 
ordinary matter is restricted to "live" on the brane. Besides this quadratic term, there is 
a dark radiation term (which may contribute with negative values for the source terms) 
as well as a cosmological constant, if a cosmological term is allowed in the bulk. Hence, 
the dynamics for our Universe is very different from that predicted by the traditional 
standard cosmological model. Due to the quadratic character of the new ordinary source 
term, its effect is in principle important during the very early Universe: The contribution 
of the quadratic term decreases very rapidly as the Universe expands. 

Many cosmological consequences may be extracted from this primordial scenario, and 
the brane cosmology became the object of intensive studies in these last times. Here, 
we will address one important point concerning the brane scenario: The primordial nu- 
cleosynthesis. In fact, the presence of a quadratic matter term and of a dark radiation 
term may have many consequences for the abundance of the primordial light elements, 
like deuterium and helium. Such problem has already been treated in references jSl E] • 
However, in reference [H] mainly the question of the influence of dark energy has been 
investigate; in reference 0, on the other hand, a balance between the contribution of 
the quadratic term and dark energy is required in order that the traditional scenario for 
nucleosynthesis is not substantially affected. These authors have kept the main features 
of primordial nucleosynthesis process, like the energy scales. In reference jH] it has been 
concluded that large negative values for the dark radiation are allowed. 

In the present work we will take two main different considerations. First, we will 
emphasize the role of the quadratic term in the source of the Einstein equations. This 
quadratic term changes drastically the temperature curve, and the characteristic energy 
scales under which the nucleosynthesis takes place becomes completely different with 
respect to the standard model. In particular, it will be shown that the observed value for 
the helium abundance implies that the nucleosynthesis is effective between 100 MeV < 
T < 0.1 MeV, in contrast with the usual values 1 MeV < T < 0.1 MeV. The fundamental 
mass scale which leads to values for the helium abundance compatible with observational 
data is of the order of M ~ 5 x 10 3 GeV. In this case, the dark radiation and cosmological 
constant may indeed be neglected at least in what concerns nucleosynthesis. In some sense, 
the analysis performed here is more restrictive than those made in references jHJ H| since 
these authors have determined a lower bound for the fundamental mass scale. However, 
we must stress that a more general analysis taking into account the contribution of all 
terms in the brane model, which asks for a substantial modification of the semi-analytical 
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method, may lead to higher values for the fundamental mass scale. 

The other aspect that distinguish the present work from the previous ones concerns 
the employment of the semi-analytical analysis for the helium synthesis which has been 
developed in reference jTUj. This semi-analytical analysis, which avoids the use of large 
numerical codes, is possible if the Fermi-Dirac statistic is ignored. Even if such statistics 
seems to be very relevant in the evaluation of the transmutation process that occurs in 
the primordial nucleosynthesis, it does not affect sensibly the final results: The semi- 
analytical method leads to a value of the helium abundance only some per cents different 
from that obtained when the Fermi-Dirac statistic is taken into account. A very important 
advantage of using such semi-analytical method concerns the fact that the main physical 
steps in the computation of the helium abundance may be emphasized. 

This article is organized as follows. In the next section, we sketch the main features of 
brane cosmology. In section 3, the semi- analytical method is exposed in its general lines. 
Its application to the brane cosmology is presented in section 4, while in section 5 the 
conclusions are presented. 



2 Brane cosmology 

The brane model to be considered here is the simplest one, which keeps the main features 
of more complicated models. This model is constructed in a five- dimensional space-time, 
where the fifth coordinate is not compact. This marks a strinking difference with respect 
to the usual Kaluza-Klein models. In this model, our usual world is constrained to a three- 
brane, where ordinary matter and fields "live". Gravity, on the other side, is defined in 
the entire five- dimensional space-time. The field equations read 

G AB = kT ab + Kg AB , (2) 

where A,B = 0,1,. ..,4, A is the cosmological constant which also "lives" in the whole 
space-time and k is the gravitational coupling in five-dimensions. The five-dimensional 
space-time is called bulk. The metric describing this configuration is 

ds 2 = n 2 (t, y)dt 2 — a 2 (t, y)^^dx l dy^ — dy 2 . (3) 

The energy-momentum tensor takes the form 

Tb = diag(p b , -p b , -p b , -p b , 0)5 (y) (4) 

where the delta function assures that ordinary matter is confined on the brane which is 
located at y — 0. 

The presence of the brane implies the need to impose match conditions through the 
extrinsic curvature. Using the field equations and considering these match conditions, we 
end up with the following equation driving the dynamics of the scale factor and matter 
on the brane: 

fd \ 2 k 4 2 A k C . . 

— = i^Pb + « - — + — > (5) 
\a J 3o o % Oq 
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the subscript indicating that all the quantities are evaluated on the brane. This subscript 
will be ignored henceforth. As usual, k denotes the curvature of the spatial section, and C 
is an integration constant generally called dark radiation, since it appears in the equation 
of motion similarly to a radiation term. Notice, however, that C may take any sign. An 
important characteristic of brane cosmology is the fact that the matter term appears 
quadratically in the equation of motion. In what follows, we neglect the curvature, fixing 
k = 0. 

The standard Friedmann-Robertson- Walker cosmology is characterized by a linear 
matter term, in opposition to equation (JSJ). However, it is possible in brane cosmology 
to generate such linear term. This can be achieved by redefining the matter term as 
Pb = p + cr, where a is a constant. Hence, © becomes 

a\ 2 k a 2 « 4 tr n A a 2 AC 

~) = ^P + T^P + + ¥ + ~ • 6 
ay 3o 18 do o a 4 

If we make the identification 8tcG = k 4 <t/18, the ordinary matter term is reproduced 
besides the quadratic one. At same time, a new cosmological term appears. The total 
cosmological term may be set equal zero by fine tuning their values. This is possible 
mainly in the context of the Randall- Sundrum model, where an anti-deSitter space-time 
is considered. 

Since we are dealing with a primordial cosmological model, the matter content of the 
universe is dominated by radiation. Hence, p = po/a 4 . In the computation that will be 
performed later, we will ignore the dark radiation, the linear term for the matter density 
and cosmological constant terms. These simplification will be justified a posteriori, since 
these terms will come out to be smaller than the quadratic matter density one for the 
energy scales relevant for the nucleosynthesis analysis. Hence, our fundamental equation 
to be analysed is 



a J M 6 900 M 6 ' 1 ; 
where we have used the fact that 

p = N—T i , (8) 

H 30 V ' 

N being the effective number of degrees of freedom of the massless particles. Since, at 
the moment of the nucleosynthesis there are photons, electron-positron pairs and three 
types of neutrinos-antineutrinos, = 43/4. The gravitational coupling parameter has 
been redefined as k 4 /36 = 1/M 6 . The parameter M defines the fundamental mass scale 
of brane cosmology. 

The use of the equation (J7|) distinguishes strongly the present approach from previous 
ones. In references jEJIH], it has been used instead of (J7J), the equation 



d\ 2 k 4 9 8%G C 



The main conclusion was that a large contribution of dark radiation is allowed if C < 0. 
But, there is one point which must be commented. In those references, the authors have 
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kept the nucleosynthesis energy scales between I MeV and 0.1 MeV. They established 
the observational limits on the dark radiation term and on the quadratic term by requiring 
that these energy scales are not substantially affected by the presence of these new terms. 
In this way, the prediction for helium abundance (as well as for deuterium and lithium) of 
the standard model, which agrees remarkably well with the observations, is not affected. 

Here, the emphasis will be different. We will use (J7J) and will leave the dynamics of 
the Universe dictates the new range of energy where the nucleosynthesis takes place. In 
this way, observational limits will be established for the term M in equation (JJJ). 

The main reason for using (J7J) comes from the following fact. The nucleosynthesis of 
helium depends crucially on the reactions converting free neutrons into proton and on the 
decay of free neutrons. The last effect continues until the energy of the Universe drops 
much below the deuterium bindind energy, e D = 2.23 MeV. Typically the decay of free 
neutrons continue until about T ~ 0.1 MeV, since the large ratio of photons to baryons 
assure that even in this energy scale there are enough photons able to dissociate the formed 
deuterium. Below this energy, deuterium is not dissociated anymore and they can combine 
to form the helium. This fact depends only on the statistic of the deuterium immersed 
in a photonic gas, depending very little on the dynamics of the Universe. However, the 
reactions converting neutrons to protons are effective until the moment the reaction rate 
is equal to the rate of the expansion of the Universe, when the number of free neutrons 
due only to this effect is frozen. 

The equality between the reaction rate TZ and the rate of the expansion of the Universe 
implies 

K ~- . (10) 

a 

Since, in the model dictated by equation (JJJ), 

(11) 




and the reaction rate for the conversion of neutrons into protons is given by 1Z = G 2 F T 5 , 
where Gp ~ 10~ 5 GeV~ 2 is the Fermi constant of the weak interactions, we find that the 
frozen temperature is given by 

' , (12) 



in opposition with the standard case where 



^3 



F n2 



GiM Pl 



(13) 



We will treat M as a free parameter, used to fit the observed helium abundance. The 
final results will justify the fact that the quadratic term in fact dominates the right hand 
side of the Einstein equations. But, before computing this, let us first describe the semi- 
analytical method to be used. 
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3 The semi-analytical computation of the primordial 
nucleosynthesis 

During the thermal history of the Universe, a quark-hadron transition occurs at a tem- 
perature higher than T ~ 100 MeV. The primordial nucleosynthesis begins after this 
phase transition. Initially, the quantity of neutrons and protons is equal. As the Universe 
expands, the temperature drops, and reactions involving the neutrons, protons, neutrinos 
and electrons convert neutrons into protons. At same time, since the neutrons are free, 
and unstable, they also decay into protons. The reactions converting neutrons into pro- 
tons ends when its rate is equal to the rate of the expansion of the Universe. However, 
the decay of free neutrons continues until the energy is low enough in order the neutrons 
to be captured forming deuterium, from which the helium is formed. The quantity of 
helium produced in this process depends essentially on the quantity of neutrons that has 
survived up to the moment they are captured to form deuterium. The detailed analysis 
of all these process is quite involved, requiring the use of numerical codes to evaluate the 
different transmutation process. However, semi-analytical expressions can be worked out 
if the Fermi-Dirac statistic is neglected. This has been done in reference ^Uj, leading to 
values for the helium abundance with an error of some percents compared with the precise 
numerical calculation. Precision is slightly lost by using this semi-analytical analysis, but 
on the other hand, the physical meaning of the computation becomes more transparent. 

We will summarize the main steps and the relevant quantities. For details, the reader 
is invited to address himself to the reference ^U]. We will consider first the standard 
cosmological model, deep in the radiative phase. The ratio of neutrons with respect to 
the total baryon number as function of the temperature is 

*< r > - ^frhm (14) 

where n n (T) and n p (T) are the numbers of neutrons and protons, respectively. The main 
equation controlling this quantity is 

^ = V(t)(i-^))-MW) - (is) 

where X pn and X np are the rates of conversion of protons into neutrons and neutrons into 
protons respectively. The main process concerned are 

X np = X(u + n — > p + e~) + A(e + + n — > p + v) + X(n — > p + v + e~) . (16) 

As an example, the first one is given by 

\( v + n^p + e -)=A dp u plp e E e {l- f e )f u , (17) 

Jo 

where A is a coupling constant, the p's denote the momenta of each particle involved in 
the process, E the energy and the /'s represent the Fermi-Dirac statistic factor. The last 
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process in equation represents the neutron decay and it is not considered in a first 
evaluation. Later, the final results will be corrected taking it into account. 

In reference [TO] the evaluation of the first two rates in equation is simplified by 
approaching the Fermi-Dirac statistics by the Maxwell-Boltzmann one. This is justified 
by the fact the temperatures concerned at the moment these process take place are smaller 
than the energies of the particles. After a lengthy evaluation of all process, we end up 
with the following expression for the neutron abundance factor: 



X{y) = X eg {y) + f dy'e*'Xl q {y') exp[K(y) - K(y')] 

j 

with the following definitions: 



18) 



XeM 

K(y) 



1 + ev ' 
4 



3 +- + 

yi yZ y 



1W4 



yO y. 



45 



V2 M n 



rAm 2 



AAr(Am) 



y 



Am 
~T~ 



(19) 
(20) 

(21) 



where Am = 1.294MeV r is the mass difference between neutrons and protons, a = 255 is 
a pure number and r is the neutron lifetime. The fraction of neutrons to baryons at the 
end of all those process, X, is obtained by making y — ► oo (T — > 0), X = X(y — > oo). 
The fact that the integration is performed until T = 0, is not important since equation 
(JIEJ) has yet saturated at y ~ 10, that is, T ~ 0.1 MeV, the temperature for which the 
neutrons were already captured forming deuterium, and the nucleosynthesis process has 
finished. 

As stated before, initially the neutron decay is neglected. To correct the final results 
due to it we must evaluate the time where the capture process forming the deuterium 
occurs. This time is obtained, in principle, by evaluating the capture process until a 
temperature of the order of the deuterium binding energy t£> ~ 2.225 Me V. However, 
lower energies (of order of E ~ 0.1 MeV) must be considered due to the fact that the 
enormous number of photons implies that deuterium dissociation continues to occur even 
when T 7 < ei). In peforming this analysis, we must take into account the evolution of the 
Universe, which in this case is reflected by the equation 



t 
_2 



inp 



3MJ, 



1/2 



(22) 



since during the radiative phase a oc 1/T 7 , T 7 being the photon temperature, which 
is approximately equal to the neutrino temperature at the relevant temperature scales 
considered in this computation. Moreover, the energy density is given by 



P 



30 



where 



N, 



eff 



UN 4/3 



iVl ~ 13 



(23) 



(24) 
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The equation (|22j) allows to establish a relation between temperature and time. In 
doing this we can convert the expression for the neutron decay in terms of the temperature. 
The final results imply that time is related with temperature by 



45 iVarii 



16vr 3 iV, 



eff 



2 / 3 M F 



T 2 

J 70 



The constant to is 
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6N, 



eff 



;^) 1/3 -i 

v 4 1 



U 



45 



(25) 



16vr 3 iV, 



eff 



T 2 
1 70 



(26) 



where T 70 is the photon temperature when the neutron capture is accomplished. The time 
of capture is obtained through some statistic considerations, and its expression is given 
by solving the equation 

2.9 x 10-V 17/6 e- L44 ^ /3+ ^ ~ 1 



where z c = e^jT. The final abundance is given by 

X f = exp {-t c /r)X 



(27) 



(28) 



where t c is the time when the neutrons are captured into deuterium. The final helium 
abundance by weight ratio is 

Y 4 = 2X f . (29) 

The complete computation yields Y4 ~ 0.247. A more precise computation, using numer- 
ical code, gives Y 4 ~ 0.241 |TT| . 



4 Nucleosynthesis in the brane cosmology 

All the computation of the helium abundance in a brane cosmology definite by equation 
(J7J) is the same as before. The main changes concerns the computation of the function 
(l20|) and the relation (|2"7j) due to the employment of equation instead of the tradi- 
tional Friedmann equation. Due to the different dynamics of the Universe, driven now by 
equation (JJJ), the function (ffllj) is now replaced by 

K'iy) = —b'\ y (1 + e~ v ) - e~ y + 6\ny + 6 / — dt\ , (30) 

I y J y t J 

where b' = a ^2 ; all other quantities remaining the same. An important feature of 
expression (p5U|) is that it does not saturate as in the standard case, due to the presence 
of the linear and logarithmic terms. Hence, the integration must be performed until 
the capture time of neutrons into deuterium, and not until T = 0. The behavior of 
the function (J30j) . compared with the corresponding function for the standard model, is 
displayed in figure (JTJ) . 

Note that the transmutation of neutrons and protons in the brane cosmology remains 
effective longer after the frozen temperature since the ratio between the expansion of the 
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. 15 



Figure 1: The behavior of the function K(y) for the standard model and for the brane 
model until the capture time. 



Universe and the reaction rate of weak interactions scales, in the brane model, with T, 
while in the standard model this ratio scales as T 3 . 

On the other hand, the capture time is also modified by the dynamics of the Universe, 
since relation (|22j) is no longer valid, being replaced by 



21 Nn 2 T 4 



T 7 30 M 3 



(31) 



This leads to the following new expression for the capture time: 



(32) 



where 



t = 37rQ Y m p \ 

\2ixae D ) 



1/3 



(33) 



with a being the fine structure constant and m p the proton mass. 

All the computations depend on the value of the parameter M. The results are 
displayed in table 1. 
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l/M 3 (MeV~ 3 ) 


X 


Y 


10 20 


0.5 


1 


10 5 


0.5 


1 


1 


0.5 


1 


10- 5 


0.5 


1 


10" 12 


0.5 


1 


10- 13 


0.5 


0.999999 


1Q -14 


0.499999 


0.999996 


7.90 x 10" 21 


0.233139 


0.250547 


7.38 x 10" 21 


0.223207 


0.230567 


1Q -4U 







Table 1: Helium abundance in terms of the mass scale value. 

From the table 1, it is possible to verify that until a value M 6 ~ 10 24 MeV 6 the value of 
the five dimensional coupling parameter is so high that it drives a so fast expansion that no 
transmutation occurs. Hence, the nucleosynthesis period finishes with the same quantity 
of neutrons it begins, that is, n p = n n . All neutrons and protons are used to produce 
helium, what explains the result Y — 1. On the other hand, for values of M 6 ~ 10 80 , the 
expansion is so slow, that no neutrons are left: The number of helium produced is zero. 
The values for the five dimensional coupling parameter such that the observational limits 
for the helium abundance are satisfied imply 5.0 x 10 3 GeV < M < 5.1 x 10 3 GeV, giving 
a very good agreement with the observed value, which is 0.23 < Y° bs < 0.25 [12J. 

Now, we can justify the approximations fixed before. First notice that the frozen 
temperature for the best values of the mass scale is now Tp ~ 75 MeV, that near the 
quark-hadron phase transition temperature. This is not in principle a problem if this phase 
transition occurs really a little before that energy. Let us consider now the equation ©. 
With those best values for M, and making the identification 



K 



2a 



8nG 



57T 



M 2 Pl 



36 M 6 ' M 6 

we have at the temperature of Tj ~ 100 MeV, the following values for each term: 
,2 



Pr 
M 6 



6 x 10" 25 MeV 



Pr 



6 x IQ- 37 MeV 



A ~ W~ 49 MeV 



(34) 



(35) 



Hence, under these conditions, the two last term in the right hand side of equation (jUJ) can 
indeed be neglected. On the other hand, the nucleosynthesis process ends when the energy 
is about Tf ~ 0.1 MeV. At this energy, the terms on the right hand side of equation © 
read 



Pr 

M 6 



6 x 10" 49 MeV 



Pr 



Ml, 



exlQ-^MeV" 2 , A~10~ 49 Me1/ 2 . (36) 



This assures that the quadratic term really dominates over the linear and the constant 
term during all the nucleosynthesis period, justifying a posteriori our first approximation. 
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Concerning the dark radiation term, we may consider that it can also be ignored during 
the nucleosynthesis period since its energy density today must not exceed too much the 
radiation density today. It scales like the ordinary radiation density; hence, its value 
at the nucleosynthesis period must be at most of the order of the linear matter term, 
much smaller than the quadratic term at least until T ~ 0.1 MeV. Notice however, that 
the constant term leads to a vacuum energy today of the order py ~ 10~ 16 GeV 4 , much 
higher than the observed value p^ 5 ~ 10~ 47 GeV 4 [13] . Even if this cosmological term 
does not spoil the nucleosynthesis process, it is not compatible with the later evolution 
of the Universe. The model must be supplemented with a mechanism to circumvent this 
problem through a dynamical process or a fine tuning. 

5 Conclusions 

In this paper we analyzed the primordial production of helium using a semi-analytical 
approach in the context of brane cosmology. The main point in this computation was to 
consider the predominant role of the quadratic matter density in brane cosmology, since 
this term dominates at very high energy scales. In doing so, we ignored the contribution 
of dark radiation. This assumption was justified a posteriori, since dark radiation and 
the ordinary linear radiation term are really a subdominant component to energy scales 
down to T ~ 0.1 MeV, if the fundamental mass scale is such that M ~ 5 x 10 3 GeV. 

The main difference with respect to previous work on the subject concerns the compu- 
tation method and the energy scales where the primordial nucleosynthesis occurs. In what 
concerns the computation method, we have opted for a semi-analytical method developed 
in reference [10J, which slightly sacrifies the precision but furnishes a more physical in- 
sight. Concerning the energy scales for the nucleosynthesis, we allow them to vary. In 
references |E1IH], the authors kept the energy scales as the same as in the standard cosmo- 
logical model, that is 1 MeV < T < 0.1 MeV. Then, they constrained the dark radiation 
term and the quadratic matter density term such that the predicted helium abundance 
remains inside the observational limit. Here, on the other hand, the energy scales for the 
primordial nucleosynthesis were kept free, and we determine them, and at same time the 
fundamental mass scale, in such a way that the observational data can be reproduced. 

The main results can be resumed as follows. If the mass scale is to small, there is no 
time for the transmutation of neutrons into protons and all baryons finish to form helium; 
no or almost no hydrogen is left. This, of course, contradicts observations. However, if 
the mass scale is too high, essentially all neutrons are transmuted into protons and no 
helium is produced, contradicting again the observations. There is an optimal value for 
the mass scale, where the observational abundance of helium can be produced. For this 
optimal value, the frozen temperature for the transmutation of neutrons into protons is 
of the order of Tp ~ 75 MeV, while the neutrons capture time inside deuterium remains 
essentially the same as in the standard scenario, i.e, T ~ 0.1 MeV. The value for the 
mass scale is M ~ 5 x 10 3 GeV. 

These results are quite reasonable. With these values, the initial hypothesis that the 
contribution of dark radiation, the linear matter density term and of the cosmological 
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constant are negligible are perfectly justified. However, there are two main difficulties 
arising from this analysis. First, the frozen temperature is dangerously near the typical 
temperature for the quark-hadron phase transition, which is not modified by the change in 
the dynamics exactly because it is a phase transition. But, there are yet many discussions 
on the exact value of the temperature at which this phase transition occurs jT3] . Second, 
the value of the cosmological constant, even if unimportant for the nucleosynthesis process, 
is too high in view of current observational limits on the vacuum energy. Some fine tuning 
is necessary in order to avoid problems in this sense. But this seems to be a common 
feature of many brane models. It is important to stress, however, that the analysis 
performed here is somehow a simplification of the problem: A more complete analysis 
should taken into account the contribution of all terms, constraining the mass scale in 
presence of the linear density term and dark energy. But this approach asks for important 
modifications of the semi-analytical method employed here. In any case, we have shown 
that a scenario where the quadratic term dominates during nucleosynthesis is possible, 
provided that the relevant energy scales are modified. 

In our point of view, the main contribution of this analysis concerns the value pre- 
dicted for the mass scale. It is somehow below other limits. In reference [S] for example, 
the authors have determined a bound such that M > 13 TeV. Our results are marginally 
consistent with this bound and, moreover, they set the fundamental mass scale of branes 
near the terrestrial accelerators energy limit, and very near the scales of other fundamen- 
tal interactions. If we had taken into account the other contributions, we could obtain a 
higher value for the mass scale. In view of a more general analysis, the result found here 
for the mass scale may be seen as lower bound. 
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